High-temperature spin relaxation process in Dy 2 Ti 2 07 probed by 47 Ti-NQR 
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We have performed nuclear quadrupole resonance (NQR) experiments on 47 Ti nuclei in Dy2Ti207 
in the temperature range 70 - 300 K in order to investigate the dynamics of 4/ electrons with strong 
Ising anisotropy. A significant change of the NQR frequency with temperature was attributed to 
the variation of the quadrupole moment of Dy 4/ electrons. A quantitative account was given 
by the mean field analysis of the quadrupole-quadrupole (Q-Q) interaction in the presence of the 
crystalline-electric-field splitting. The magnitude and the temperature dependence of the nuclear 
spin-lattice relaxation rate was analyzed, including both the spin-spin and the Q-Q interactions. 
The results indicate that these two types of interaction contribute almost equally to the fluctuation 
of Dy magnetic moments. 
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I. INTRODUCTION 

Observation of spin relaxation in magnetic materials 
provides valuable information on the microscopic nature 
of the interactions between spins. When a system of a 
large spin quantum number has strong Ising anisotropy, 
the mutual spin-flip process is substantially depressed. 
The relaxation then becomes extremely slow and may 
eventually lead to spin freezing in real materials with 
competing interactions. A typical example is Dy 2 Ti 2 07, 
in which Dy 3+ ions with large angular momentum of 
J=15/2 (4/ 9 configuration) form a pyrochlore lattice, a 
network of corner-shared tetrahedra. 

InDy 2 Ti 2 7 (Ref.H), as well as in other pyrochlore ox- 
ides such as Ho 2 Ti 2 07 (Ref.[2|) and Ho 2 Sn 2 C>7 (Ref. y), a 
strong crystalline electric field (CEF) forces the magnetic 
moments of rare earth ions to lie at low temperatures 
along the symmetry axis of CEF, i.e. the [111] or its 
equivalent directions connecting the rare earth site and 
the center of the tetrahedron it belongs to. This strong 
Ising anisotropy combined with the effective ferromag- 
netic nearest neighbor interaction stabilizes for each bond 
such magnetic configurations that one moment points in 
and the other points out of the tetrahedron. Obviously it 
is not possible for all bonds to satisfy this condition. As a 
compromise, the ground state is selected by the rule that 
two moments point in and other two moments point out 
of every tetrahedron, which is called the "two-in, two- 
out" or the "ice rule". These states are macroscopically 
degenerate as long as only the nearest neighbor interac- 
tion is considered. This state has been named "spin-ice" 
after the Ih phase of water ice, which exhibits similar 
degeneracy with respect to the proton configuration. 

The process of spin relaxation in such a system, es- 
pecially in the spin-ice state, is an interesting problem. 
Although the degeneracy of the spin-ice state may be 
eventually lifted by the long range dipolar interaction, 
this does not occur in real materials within an observ- 
able time scale. Thus the true ground state may never 



be reached. Recently, the relaxation in Dy 2 Ti 2 07 and 
Ho 2 Ti 2 07 have been actively investigated by various ex- 
perimental means, including the neutron spin echo ex- 
periments in Ho 2 Ti 2 07 (Ref. 0), the muon spin rotation 
measurements in Dy 2 Ti 2 7 (Ref. d) and the observation 
of magneto-caloric effects in Dy 2 Ti 2 07 (Ref. 0). How- 
ever, the time scale of magnetic relaxation in Dy 2 Ti 2 07 
and Ho 2 Ti 2 07 varies with temperature so widely that no 
single experimental technique can cover the whole range. 

In this paper, we report the results of nuclear mag- 
netic resonance (NMR) and nuclear quadrupole reso- 
nance (NQR) experiments on 47 Ti nuclei in Dy 2 Ti 2 07. 
From the analysis of the NMR spectra in high magnetic 
field and at low temperatures, where the configurations of 
Dy spins are precisely known, we conclude that the mag- 
netic hyperfine field at Ti nuclei is entirely due to dipo- 
lar interaction with Dy moments. The zero-field NQR 
measurements have been performed in the temperature 
range 70 - 300 K. Substantial temperature variation of 
the NQR frequency was attributed to the change of the 
quadrupole moment of Dy 4/ electrons. A quantitative 
account of the data was obtained by considering both 
the CEF splitting and the quadrupole-quadrupole (Q-Q) 
interaction in a self-consistent way. The fluctuation rate 
(1/t) of the Dy magnetic moments extracted from the 
nuclear spin-lattice relaxation rate shows an activated 
temperature dependence. The magnitude and the In- 
dependence of 1/t was successfully explained by consid- 
ering both the spin-spin and the Q-Q interactions, which 
have comparable contributions. 



II. EXPERIMENT 

For the 47 Ti-NQR/NMR experiments we have grown 
an isotope-enriched single crystal of Dy 2 Ti 2 07 using a 
floating-zone image furnace, starting from the sintered 
polycrystal containing 96% 47 Ti while the natural abun- 
dance of 47 Ti is 7.5%. The purpose of this enrichment is 
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FIG. 1: (Color online) (a) The crystal structure ol Dy2Ti2C>7 
viewed perpendicular to the [111] direction. A magnetic field 
applied along [111] generates two spectroscopically inequiv- 
alent Ti sites; Ti 1 on the triangular plane and Ti 2 on the 
Kagome plane of the Ti pyrochlore sublattice. (b) 47 Ti-NMR 
field-swept spectra at 9.7 MHz for the two field orientations. 
The "one- in, three-out" configuration is selected for H || [111], 
whereas the "two-in, two-out" configuration remains stable 
for H || [001]. The arrows indicate the calculated resonance 
positions as described in the text. For H | [001], the two 
resonance lines in low field region were not observed probably 
because of too short T2. 



to enhance the signal intensity and to avoid overlap of res- 
onance lines from 47 Ti and 49 Ti, which have nearly the 
same gyromagnetic ratio ( 47 - 49 -f/2n ~ 2.400 MHz/T). 
Since the demagnetization field of Dy 2 Ti 2 07 is large in 
high fields and at low temperatures, of the order of one 
tesla, the single crystal was cut into a rod extending along 
the [110] direction and the NMR coil was wound only 
around the central half of the rod where demagnetiza- 
tion field should be homogeneous. 

The standard spin-echo technique was used to obtain 
the NMR/NQR spectra. The nuclear spin-lattice relax- 
ation rate (1/Ti) was measured by NQR using the satu- 
ration recovery method. To determine 1 /Ti , the intensity 
of the spin-echo signal as a function of the time after the 
saturation pulse was fit to the theoretical formula ap- 
propriate for the specific resonance line being observed^. 
The spin-spin relaxation rate 1/T 2 was determined by 
fitting the spin-echo intensity as a function of the time 
interval (r) between the tt/2 and the tt pulses to an ex- 
ponential function. 

Since the NQR frequencies are relatively low (4 and 
9 MHz), the long decay time of the electric ringing of the 
LC resonance circuit after the rf-pulses prevents obser- 



TABLE I: Comparisons between the experimental and calcu- 
lated resonance positions for 47 Ti NMR at 9.7 MHz. The cal- 
culated results are obtained by exact diagonalization of 7i nU c 
[Eq. (JTJ)] with or without considering the hyperfine field from 
the Dy moments ff dip [Eq. @]. 
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vation of the spin-echo signal for short r. This raises a 
serious problem at low temperatures where T 2 becomes 
short. In order to overcome this difficulty, we have devel- 
oped a silicon-diode-based Q-switch. This device quickly 
dumps the quality factor (Q) of the resonance circuit just 
after the irradiation of rf-pulses, allowing us to set r as 
short as to 8 /Lis at 9 MHz. 

Still the NQR experiments are not possible in a cer- 
tain temperature range because of too short T 2 due to 
the slow dynamics of Ising spins 8 . The measurements 
reported in this paper are limited to high temperatures 
(> 70 K), where the spin dynamics is fast enough, or at 
low temperatures (< 4 K) and in high fields (> 3 T), 
where magnetic moments are completely aligned to cer- 
tain directions by magnetic fields. 

III. RESULTS AND DISCUSSION 

A. High field NMR 

The crystal structure of Dy 2 Ti 2 0y shown in Fig. [IJa) 
contains two pyrochlore sublattices, one formed by Ti 
and the other formed by Dy. A pyrochlore lattice can be 
viewed as the stack of alternating triangular and Kagome 
planes along [111]. When a magnetic field applied is 
along [111], The Ti sites on the triangular planes [Ti 1 
in Fig. Ufa)] and those on the Kagome planes (Ti 2 ) be- 
come inequivalent, generating distinct NMR lines. Note 
that Ti 1 (Ti 2 ) is on the Kagome (triangular) planes of the 
Dy pyrochlore sublattice. When the field is applied along 
[100], all Ti sites are equivalent. Figure HJb) shows the 
NMR spectra obtained at the fixed frequency of 9.7 MHz 
for the two field orientations, along [111] and [100], at low 
enough temperatures. The magnetization measurements 
have established that at these fields and temperatures, 
each Dy moment of 10/xb is completely aligned to a spe- 
cific directio n 9 ' 10 ! 11 ! 12 as shown in the inset of Fig. [IJb). 
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The "one-in, three-out" or the "three-in, one-out" con- 
figuration is selected for H | [111], while the "two-in, 
two-out" configuration remains stable for H || [100]. In 
the following, we demonstrate that the NMR spectra can 
be well reproduced by considering the dipolar field from 
Dy moments in these configurations. 

The NMR frequencies of 47 Ti are determined by the 
following nuclear hamiltonian, which consists of the Zee- 
man interaction with the effective magnetic field and the 
nuclear quadrupole interaction with the electric field gra- 
dient (EFG)±2, 
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Here J is the nuclear spin with 7=5/2, Q (= 3.02 x 
m 2-j } g nuc i ear quadrupole moment and V a ° (= 
d 2 V/dr a drp) is the EFG tensor, where V is the electro- 
static potential. The repeated indices of vectors and ten- 
sors should be summed over the Cartesian coordinates 
(Einstein's convention). 

The effective field H e s is the vector sum of the ex- 
ternal field H and the hyperfine field. We assume that 
the hyperfine field is given by the dipolar field from Dy 
moments. 
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where hb is the Bohr magneton, gj (= 4/3) is the Lande's 
g factor, Id is the demagnetization factor, and Va is the 
atomic volume per Dy. The direction of the saturated 
Dy moments J (= 15/2) are fixed along the local Ising 
axes in the low-temperature limit of our experiment. The 
summation over the spins j should be cut off within a 
sphere or performed by the Ewald's method. We take 
7d=0.45, which is close to the value 1/2 for a infinitely 
long cylinder in horizontal fields. For H || [111], the 
calculated dipolar field is (-0.30, -0.30, -0.30) T at Ti 1 
sites and (-0.24, -0.24, -0.66) T at Ti 2 sites. For H || 
[001] all Ti sites have the same dipolar field (0.18, 0.18, 
-0.72) T. 

The EFG tensors can be determined from the NQR 
measurements at zero magnetic field. Since the Ti sites 
have three- fold rotation symmetry around [111], the EFG 
is axially symmetric, V xx = V yy — —V zz /2, where z \\ 
[111]. When H c h = 0, Ttnuc can be easily diagonalized to 
yield two NQR frequencies, v Q [= 3eQ\V zz \/2I(2I - l)h] 
and 2isq, corresponding to the transitions I z = ±1/2 
±3/2 and I z = ±3/2 ±5/2, respectively. 

Now we can diagonalize 7i nuc numerically to obtain the 
NMR frequencies, which are the difference of the eigen- 
values of Ti nuc for various transitions, as a function of 
external field. In order to compare the calculated results 



with the experimental NMR spectra obtained by sweep- 
ing the field at a constant frequency, this function should 
be inverted by Newton's method. The experimental and 
calculated results for the NMR line positions are com- 
pared in Fig.QJb) and Table HI Here we used the value of 
vq = 4.39 MHz to get the best agreement. This value is 
very close to the observed NQR frequency at low temper- 
atures as we describe below (Fig. [3]). The results without 
the hyperfine field are also shown in Table HI The calcu- 
lated results including the dipolar field agrees extremely 
well with the experimental data, demonstrating that the 
hyperfine coupling in Dy2Ti2 07 is entirely due to the 
dipolar field. The precise knowledge of the hyperfine in- 
teraction enables us to make quantitative analysis of the 
NQR relaxation data we discuss later. 



B. Zero field NQR 

We now focus on the NQR results in zero field and at 
high temperatures. Figure[5]shows NQR spectra at 240 K 
and at 150 K. Two resonance lines are observed at vq 
and 2vq as expected. However, the value of vq strongly 
depends on temperatures, which is rather unusual. The 
temperature dependence of the peak frequency and the 
width of the NQR lines are shown in Fig [3] The value 
of isq changes as much as 5% from 300 K to 70 K. This 
is an order of magnitude larger than what is expected 
from simple thermal contraction^. The smooth change 
of vq and the nearly constant line width also rule out any 
sudden deformation of the crystal structure. A natural 
explanation would be that vq is influenced by the change 
of charge density distribution of Dy 4/ electrons due to 
variation of the population among the CEF split 4/ lev- 
els. For example, Tou et al. reported that the tempera- 
ture dependence of EFG at Sb nuclei in PrOs4Sbi2 can 
be accounted for by the thermal average of the hexade- 
capole moments of Pr 2+ (4/ 2 ) ions which is compatible 
with the cubic point symmetry of the Pr sitesi^. 

The CEF Hamiltonian Hcef of single Dy ion in 
the Z?3d symmetry is given by the CEF parame- 
ters and the normalized spherical harmonics C^[= 



150 K 
240 K 



±l/2<->±3/2 



±3/2 <-> ±5/2 



4 4.5 5 5.5 6 6.5 7 7.5 

Frequency (MHz) 



FIG. 2: (Color online) 47 Ti-NQR spectra at 150 K and 240 K. 
The NQR frequencies depend strongly on temperature. 
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mal average (C*). Because of the 1)3^ symmetry of the 
Dy ions (C*)=0 for g 5/ 0. Then to the lowest order, 
isq = a + a2(CQ)+a4(C^)+ae{CQ). Although the micro- 
scopic mechanism of the coupling is not understood, 
it should include the effect of polarization of the sur- 
rounding medium. In the following we consider only the 
E; Wave function /, C?o 
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FIG. 3: (Color online) The peak frequency (lower panel) 
and the full width at half maximum (FWHM, upper panel) 
of the 47 Ti NQR resonance lines for the transitions I z — 
±1/2 ^ ±3/2 (vq) and I x = ±3/2 ^ ±5/2 (2u Q ). The 
dotted and double dashed lines in the lower panel show the 
calculated results considering the temperature dependence of 
the quadrupole moment (Co) without and with the electric 
quadrupole-quadrupole interaction, respectively. 



FIG. 4: (Color online) CEF level scheme of Dy 2 Ti 2 7 inferred 
from the CEF parameters of Ho2Ti2C>7 fRef. KLq) . Arrows rep- 
resent the major process for the spin-flip relaxations between 
the ground-state doublet. 



quadrupole moment and neglect the higher order multi- 
poles, 



v /47r/(2fc ± l)Yq] as operators. 
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The values of B* 



for Ho2Ti207 have been determined by 
Rosenkranz and co-workers^ by neutron scattering ex- 
periments. Among the family of lanthanoid titanate py- 
rochlorcs, should change only modestly We estimate 
B^ in Dy2Ti207 by assuming the change of (r k ) for the 

4/ shells as (r^n y /(r^) H o= 1-05, 1.09, 1.12 (for fc=2, 4, 
6), after Refs.[r7l and Ho. In order to calculate the matrix 
elements of C*, the \ J Z ) states have to be expanded to 
the single-electron orbitals with the Clebsch-Gordan co- 
efficients. The CEF level scheme for Dy2Ti20y obtained 
by diagonalizing 7Ycef for the case of e H 15 / 2 multiplet is 
shown in Fig. [4j Our results are in good agreement with 
those described in Ref. [HI- Note that the CEF ground 
states are almost the pure \ J Z = ±15/2) doublet, leading 
to the strong Ising anisotropy. 

Since the charge distribution of 4/ electrons is gener- 
ally represented by the electric multipole moments C„, 
it should be possible to expand Vq in terms of their ther- 



= aa + 02 (Co)- 



(4) 



If we do not consider the quadrupole-quadrupole (Q-Q) 
interaction between neighboring Dy sites, (Cq) is simply 
given by the Boltzman average of the expectation values 
of Cq for the CEF eigen states, 



(C 2 ) 



(5) 



where/? = (ksT) -1 , and Z = ^2 i exp(—(3E i ). The exper- 
imental data of i/q is fit to Eqs. ^ and ((3J) with the two 
adjustable parameters clq and a-i as shown by the dotted 
line in Fig. [3] (labeled as "CEF calc. w/o Q-Q"). Al- 
though the experimental data are reproduced quite well 
above 120 K, clear deviation develops at lower temper- 
atures. Modest change of the CEF parameters does not 
lead to notable improvement. In the following we demon- 
strate that substantial improvement can be achieved by 
considering the electric quadrupole-quadrupole (Q-Q) in- 
teraction in the mean field approximation. 

The traceless quadrupole moment Q a p in the Carte- 
sian coordinate is defined as^, 
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In this definition, effects of shielding is already taken into 
account by the Sternheimer shielding factor cr 2 = 0.527 
and (r 2 ) = 0.849 a. u3. From the multipole expansion of 
the Coulomb potentials^, the Q-Q interaction between 
two Dy ions is expressed as 
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It should be noted that in the formulas of the Cartesian 
coordinates, the operators defined on the local quantiza- 
tion axes have to be converted as RQR T and RJ, where 
R is the rotation matrix. 

The single site Hamiltonian Ti^ 1 ' including both CEF 
and the Q-Q interaction in the mean-field approximation 
can be written as follows. 



CEF 



(0 
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167T60 



(gg) (gS) 

dr a drpdr 7 dr s \r tj J (C%) (C$) 



= Uc,,'-'+(96 K)x [C 2 ^\CI)-\{CI) 2 



(8) 



The thermal average value of (Cq) is calculated again 
from Eq. |5| but now Ei and (Cg)i are the energy and 
the expectation value for the eigenstates of the mean field 
hamiltonian, Eq. ([8]), which itself contains (Cq) as a pa- 
rameter. The temperature dependence of (Cq) can then 
be determined self consistently by iteration. Using this 
result, the experimental data of vq is fit to Eq. as 
shown by the double dashed line in Fig. [3] (labeled as 
"CEF calc. w Q-Q"). We notice that the calculated re- 
sult agrees with the experimental data much better now, 
in particular below 100 K. Although the Q-Q interac- 
tion is necessary to consider the high-temperature excita- 
tions, the spin-spin interaction is much more important in 
the low-temperature physics since the ground-state (GS) 
Kramers doublet has no freedom of multipole moment. 
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FIG. 5: (Color online) Temperature dependences of the spin- 
lattice relaxation rate 1/Ti and the spin-spin relaxation rates 
I/T2 measured for the I z = ±1/2 <-» ±3/2 (uq) and I z = 
±3/2 <-> ±5/2 (2fQ) transitions. The 1/Ti was measured at 
2uq. Three straight lines indicates the T -3 dependence. One 
is the fit to the data of 1/Ti and others are multiplied by 4.6 
and 7.6 respectively (see text). 



C. Spin dynamics 

Let us now discuss the results of the nuclear spin-lattice 
relaxation rate 1 /T\ and the nuclear spin-spin relaxation 
rate I/T2 measured by NQR. They provide direct infor- 
mation on the fluctuations of Dy moments at the NQR 
frequencies. Figure [5] displays the temperature depen- 
dence of 1/Ti and 1/T 2 . For the whole temperature 
range, l/Tj obeys the 1/T 3 behavior remarkably well. 
Similar temperature dependence is observed for 1 /T 2 be- 
low 200 K. For both 1/Ti and 1/T 2 , the relaxation pro- 
cess is considered to be magnetic because relaxation rate 
by phononic possesses should increase with increasing 
temperature 2 ^. 

The spin-spin relaxation rate is the sum of two con- 
tributions, 1/T 2 = (l/T 2 ) Dy + (1/T 2 ) nuc . The first term 
is due to fluctuations of Dy moments and the second 
term comes from the nuclear spin-spin coupling, When 
the fluctuations of Dy moments are much faster than the 
NQR frequencies, (l/T 2 )D y for the transition I z = m «-> 
to + 1 is related to 1 jT\ as 



=i{I(/+l)-m(TO+l)-l} + -L, (9) 



This relation can be derived by extending the Walstedt's 
method for the transition I z = 1/2 -1/2 (Ref. HI). 
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Here T| (T-j 1 ) is the relaxation time when the nuclear 
spin is quantized parallel (perpendicular) to z, the prin- 
cipal axis of EFG. Note that T-J' is equal to the actual re- 
laxation time measured by NQR experiment. The ratio 
T-j'/Tf 1 - is determined by the anisotropy of the hyperfine 
couplings between Dy moments and a Ti nucleus. 
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1A\ 



0.6. 



(10) 



Here, A\\ = 0.0067 T//x B (A ± = 0.019 T//i B ) for the 
component of the dipolar field from a single nearest Dy 
moment parallel (perpendicular) to z. Then, the ra- 
tio (T 2 -1 ) D y/Tf 1 should be 4.6 for m = 1/2 and 7.6 
for m = 3/2. In fact, the experimental data in Fig. [5] 
agree with these relation below 200 K, indicating that 
1/T 2 is dominated by the spin-lattice relaxation pro- 
cess (T^ )ny The deviation at higher temperatures is 
likely to be due to the temperature independent term 
(1/T 2 ) nuc . 

The rapid increase of 1/Ti with decreasing temper- 
ature indicates slowing down of the fluctuations of Dy 
moments. Basically this can be understood from the 
fact that the direct spin flip process is inhibited among 
the CEF GS doublet, which is almost purely the \J Z = 
±15/2) states. Thus the spin relaxation process at high 
temperatures should be driven by the transitions between 
the GS and the excited states. If the effects of the second 
and higher excited states are neglected, the probability 
of such a process should follow the Arrhenius-type ac- 
tivated temperature dependence exp(— E^/k^T), where 
E a is the energy difference between the GS and the first 
CEF excited states. In several experiments to date, the 
spin relaxation has been fitted to the Arrhenius function: 
the muon spin rotation and the the ac-susceptibility mea- 
surements (E a — 210 K, Refs. [B| and [If), and the nuclear 
forward scattering experiment (E a — 270 K, Ref. l22h . 
Although the 1/Ti data in Fig. [5] apparently does not fit 
to the Arrhenius function, we demonstrate in the follow- 
ing that this can be understood from the temperature 
dependences of both the amplitude of moment fluctua- 
tions and the fluctuation time r, and r does follow the 
Arrhenius function. 

When the fluctuations of Dy moments are much faster 
than the NQR frequencies with short enough correlation 
length, 1/Ti is given by the Fermi's golden rule, 



= ^(*V±5jMb) {Jt)r. 



Ti h? 



(11) 



Here we consider fluctuations of only the z component 
of Dy moments because of the strong Ising character of 
the GS. The temporal correlation of J z is assumed to 
be described by an exponential function exp(— t/r) and 
z' (=6) is the number of nearest neighbor Dy moments. 
The amplitude of the moment fluctuations (J 2 ) is calcu- 
lated from the CEF hamiltonian and shown in the inset 
of Fig. [5] The fluctuation time r can be estimated from 
Eq. (fTTj) by considering (J 2 ) and the experimental data 



of 1/Ti, as shown in the main panel of Fig. [6] for the 
temperature range 100 - 300 K. In Fig. [BJ t is found to 
accommodate to the simple Arrhenius function with E a 
of 409 ± 10 K, which is close to the CEF first excitation 
level of 385 K. We also show the temperature dependence 
of t below 100 K obtained by Sutter et al. from the nu- 
clear forward scattering 22 . The two sets of experimental 
results are consistent each other within a factor of three. 

The spin flip process relevant to the nuclear relaxation 
is illustrated by the arrows in Fig. First, a Dy mo- 
ment is excited from the GS to one of the the excited 
doublet. It then makes the transition to another state 
within the same excited doublet. This process is caused, 
for example, by the coupling J+ \ which leaves the 
spin state at site 1 in the GS unchanged and flips the 
spin at site 2 in the excited doublet. The latter process 
is allowed for the first excited doublet because it contains 
the | ± 1/2) components. Finally the Dy moment in the 
excited state relaxes back to the GS but final state must 
be different from from the initial state. This effectively 
accomplishes the spin flip in the GS: J z = 15/2 «-> —15/2. 
The transition within the excited doublet is required be- 
cause otherwise the final state always coincides with the 
initial state. Below the room temperature, the popula- 
tion of excited states is so small that the fluctuation rate 
is exclusively determined by the first step. 

The spin flip rate 1/2t is then given as the sum of 
the transition probability W from the GS to the excited 
states i when Ei 3> fceT. 



1 

27 



(12) 



The factor of 1/2 is necessary because a half of excited 
spins return to the original ground state without flipping. 
The transition probability Wgs— ►« is given by the golden 
rule for the mutual transition process. 



W G s^i ~ z — 
n 



exp(-f3E l ) 1 

Z irhWi^cs ' 



(13) 



where the spin at site 2 makes transition from the GS 
to an excited state i while the spin at site 1 does the re- 
verse and z{= 6) denotes the number of nearest-neighbor 
bonds. We consider the uncertainty principle that the 
energy width of the final state is given by the life time of 
the excited spin at site 2, which is Wi^os- From Eq. (fl3|) 
and the detailed balance relation, 

W GS ^ = Wi^cs exp(-0Ei), (14) 

we obtain the expression for W G S—>i and 1/t, 

W G s^^^\(GS,i\H^\i,GS)\ e -^t^, (1 
n ^JZ 



(15) 



|(GS,i|W 12 |i,GS) 



E 



expj-pEj) 

Vz ■ 



(16) 
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FIG. 6: (Color online) Spin relaxation time r as a function 
of temperature, extracted from the 47 Ti-NQR relaxation rate 
1/Ti using Eq. with the amplitude of the moment fluc- 
tuation (J 2 ) shown in the inset. The results of the nuclear 
forward scattering 2 ^ is shown for comparisons. The Arrhenius 
fit gives the activation energy _E a of 409 ± 10 K. The calcu- 
lated results of r based on the CEF scheme is shown for three 
cases, taking account of only the spin-spin interaction, only 
the quadrupole-quadrupole interaction, or both interactions 
together. The results for the first two cases coincide almost 
completely. 

It is rioted that Eq. (fl"75)) agrees with the Arrhenius be- 
havior if the the second and the higher CEF levels are 
ignored. 

Finally, we are at the stage to calculate 1/r in Eq. (JTHJ). 
The interaction between two moments 7i 12 , which is es- 
sential to relax the spins, is considered only between the 
nearest neighbors and contains of two terms: The spin- 
spin interaction Hgp in _ spin and the Q-Q interaction Hq.q 
given in Eq. l[7]). The spin-spin interaction consists of the 
antifcrromagnetic exchange (J cx /3 = —1.24 K) and fer- 
romagnetic dipolar interactions according to the Monte 
Carlo study of the "dipolar spin-ice model" 23 i 24 , 

-1/12 

'^spin-spin 

= + mM&^_ / j_\ | J{)jf) ^ 

I J 2 4tt dr a dr (i \r 12 J J p 

(17) 

The calculated results of r are compared with the ex- 



perimental results in Fig. [6] In the calculated results 
in, the contributions from the the spin-spin interaction 
and the Q-Q interaction turn out to be nearly the same. 
When both interactions are considered, the calculated 
result of r ("Q-Q+spin-spin" in Fig. [6]), agrees with the 
experiment reasonably well. However, there is still about 
a factor four difference. Possible origin for this may come 
from effects we have neglected such as the interactions be- 
tween second nearest and further neighbors, interactions 
between octopoles or higher order multipoles, and short 
range spin correlation. However, any of those effects is 
difficult to estimate precisely at this moment. 

In our formulation thus far, the direct-transition pro- 
cess between the GS is ignored just because such prob- 
ability is negligibly small compared with the process we 
discussed in the high-temperature region. Nevertheless, 
the probability is not strictly zero owing to the existence 
of ±9/2 and ±3/2 components in the GS (Fig. gj). We 
consider that the direct spin flip process has to be respon- 
sible for the excitations in the spin-ice state of Dy2Ti20y 
below 1 K, which is still not uncovered experimentally. 



IV. CONCLUSION 

We have investigated the Dy Ising-spin relaxation pro- 
cess in a range of 70 - 300 K using 47 Ti-NMR/NQR in 
Dy2Ti 2 07. The hyperfme coupling, CEF level scheme, 
and interaction Hamiltonians, were deduced from the 
low-temperature NMR spectra, experimental CEF of 
H02T12O7 (Ref. [HI), and the temperature dependence 
of the NQR frequency, respectively. Those information 
enabled us to calculate the absolute values of the Dy spin- 
flip rate and NQR relaxation rate. The absolute value of 
the NQR relaxation rate was analyzed in the CEF level 
scheme with two types of mutual spin- flip Hamiltonians: 
the spin-spin and Q-Q interaction. The Q-Q interaction 
is also necessary to account for the feature of the mod- 
ification in the NQR frequency. Our experimental and 
calculated results reveals the spin-flip process and the im- 
portance of the Q-Q interaction in the high-temperature 
nature of the Ising-spins in Dy 2 Ti 2 07. 
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